CTmatrix =    [0.4 0.1;
               0.2 0.5];
tol = 1;
x0 = [8;9];
pData = [2.4;3.0];
maxIter = 100;
[x res] = reconART(pData,CTmatrix,tol,x0,maxIter);